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1.0  INTRODUCTION 


A  major  component  of  the  Digital  O  blique  R_emote 
Ionospheric  Sensing  program  (DORIS)  is  the  development  of  an 
automatic  oblique  ionogram  scaling  algorithm.  The  nature  of  the 
variations  that  have  been  observed  in  oblique  ionograms  collected  to 
date  has  required  the  use  of  adaptive  signal  processing  and  pattern 
recognition  techniques  in  the  development  of  this  oblique  ionogram 
scaling  algorithm.  An  extensive  data  base  of  digital  oblique 
ionograms  will  be  required  to  develop,  test,  and  verify  the 
performance  of  the  scaling  algorithm.  The  continuing  efforts  in  data 
collection  and  the  associated  parallel  effort  in  the  development  of  the 
oblique  ionogram  scaling  algorithm  is  discussed  in  this  report. 
Results  obtained  to  date  have  lead  to  a  modification  of  the  initially 
proposed  oblique  ionogram  scaling  algorithm.  These  modifications, 
discussed  in  detail  in  the  following  section,  have  produced  a  more 
robust  as  well  as  a  more  computationally  efficient  oblique  ionogram 
scaling  algorithm. 
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2.0  REFINEMENT  OF  OBLIQUE  IONOGRAM  SCALING  ALGORITHM 


The  operation  of  the  endpoint  vertical  ionogram  based 
oblique  ionogram  scaling  algorithm  presented  in  the  Scientific  Report 
No.  13,  GL-TR-89-01 84,  Oblique  Ionogram  Automatic  Scaling 

Algorithm  -  A  Status  Report,  October  1988,  can  be  represented  by 
the  block  diagram  presented  in  Figure  2.1.  The  noise  suppression/ 
trace  enhancement  filter  operation  of  this  proposed  algorithm  was  a 
Wiener  filter  that  required  knowledge  of  the  power  spectral 
densities  of  the  residual  noise/interference  signal,  Pn(cot,cof),  and  the 
leading  edge  of  the  oblique  ionogram,  Pn(cot,cof).  In  addition  the  two- 
dimensional  Fourier  transform  of  the  oblique  echo-trace  pulse  width 
generation  filter,  Hs(cot,(Of),  was  needed  to  produce  the  Wiener  filter 
H(cot,fdf)  given  as: 


Examination  of  the  behavior  of  this  filter  on  oblique 
ionogram  data  for  paths  of  1500  and  2500  km  indicated  that  the 
filter  is  very  sensitive  to  Hs(wt,<of),  that  represents  both  the  oblique 
echo-trace  pulse  width  and  pulse  shape.  Since  we  can  only 
approximate  the  oblique  echo  pulse  shape,  the  oblique  echo-trace 
pulse  width  generation  filter,  Hs(oo  t,cof),.was  removed  from  the 
Wiener  filter,  producing  a  modified  Wiener  filter  Hm(oat,cof)  that 
produces  an  optimal  estimate  of  the  entire  oblique  echo-trace  from 
the  preprocessed  oblique  ionogram.  Hm(o)t,(Of)  is: 


HrrK,CDf)  = 


PicK’wf) 


Pli«t’wf)  +  Pn(«r«f) 
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With  this  modification,  the  detection  of  the  leading  edge  of  the 
oblique  echo  traces  will  take  place  as  part  of  the  trace  identification 
algorithm.  The  modified  endpoint  vertical  ionogram  based  oblique 
ionogram  scaling  algorithm  can  be  represented  by  the  block  diagram 
presented  in  Figure  2.2. 
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3.0  EVALUATION  OF  OBLIQUE  TRACE  ENHANCEMEN  1  ALGORITHM 


The  performance  of  the  modified  oblique  ionogram  trace 
enhancement  algorithm  was  evaluated  using  two  sets  of  ionograms. 
The  first  data  set  was  obtained  during  days  90-98  of  1988.  The 
oblique  path  was  approximately  1800  km  long  with  the  transmitting 
Digisonde  256  in  Goose  Bay,  Canada  and  the  corresponding  receiving 
system  at  Millstone  Hill  Observatory  in  Westford,  MA.  The  Goose  Bay 
station  used  a  horizontally  radiating  log  periodic  antenna,  while  the 
Millstone  Hill  station  used  the  standard  Digisonde  seven  loop  antenna 
array.  Synchronization  of  the  two  stations,  which  is  required  for 
oblique  ionograms,  was  performed  manually  in  this  experiment.  The 
vertical  ionograms  were  produced  over  a  frequency  range  of  1  to  10 
MHz,  using  100  kHz  sounding  frequency  increments,  while  the 
oblique  ionograms  were  scanned  over  a  frequency  range  from  4  to 
18  MHz,  using  50  kHz  frequency  increments.  The  second  data  set 
was  obtained  during  day  301  of  1988.  The  oblique  path  was 
approximately  2500  km  long  with  the  transmitting  Digisonde  256  in 
Goose  Bay,  Canada  and  the  corresponding  receiving  system  at 
Wallops  Island,  VA.  The  vertical  ionograms  were  also  produced  over 
a  frequency  range  of  1  to  10  MHz,  using  100  kHz  sounding  frequency 
increments,  while  these  oblique  ionograms  were  scanned  over  a 
frequency  range  from  5  to  29  MHz,  using  25  kHz  frequency 
increments. 


The  2500  km  oblique  data  is  seen  in  Figures  3.1  -  3.7. 
The  original  25  kHz  oblique  ionogram  and  the  resulting  100  kHz 
ionogram,  produced  via  the  frequency  redundancy  technique  are 
seen  in  Figure  3.1.  Figure  3.2  is  an  enlarged  view  of  the  100  kHz 
oblique  ionogram.  Figure  3.3.  shows  the  resulting  output  of  the 
modified  Wiener  trace  enhancement  and  noise  suppression  filter. 
Figure  3.4  shows  the  region  of  interest  calculated  from  the  F  layer  of 
the  corresponding  vertical  ionogram.  The  synthesized  oblique  F  trace 
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Figure  3.1  Original  25  kHz  oblique  ionogram  (upper  panel)  and  resulting  100  kHz 
ionogram  produced  via  the  frequency  redundancy  technique.  Goose  Bay,  Labrador  to  Wallops 
Island,  Virginia;  October  27,  1988,  1201UT. 


Figure  3.2  Enlarged  view  of  the  100  kHz  oblique  ionogram  shown  in  Figure  2.1.  Goose 
Bay,  Labrador  to  Wallops  Island,  Virginia;  October  27,  1988,  1201UT. 
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Figure  3.3  Resulting  100  kHz  oblique  ionogram  after  processing  through  the  modified 
Wiener  trace  enhancement  algorithm  and  the  noise  suppression  filter.  Goose  Bay,  Labrador 
to  Wallops  Island,  Virginia;  October  27,  1988;  1201UT. 


Figure  3.4  Region  of  interest  calculated  from  the  F  layer  of  the  corresponding  vertical 
ionogram.  Synthesized  oblique  IF  ionogram  is  shown  in  blue.  Goose  Bay,  Labrador  to 
Wallops  Island,  Virginia;  October  27,  1988,  1201UT. 
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is  shown  in  blue  in  the  center  of  the  white  region  of  interest.  At 
present  the  algorithm  will  accept  any  signals  in  this  region,  defined 
as  all  points  within  six  range  bins  and/or  sounding  frequencies  of  the 
synthesized  oblique  F  trace,  as  potential  oblique  F  trace  data.  At  a 
later  stage  of  development  additional  information  derived  from  an 
ionospheric  model  and  information  based  on  the  total  previous 
experience  of  the  oblique  scaling  algorithm,  i.e.  the  results  of  all 
previously  scaled  oblique  ionograms,  will  be  used  to  determine  the 
region  of  interest.  Figure  3.5  shows  the  region  of  interest  overlaid  on 
the  filtered  oblique  ionogram  data.  Note  that  although  all  of  the  high 
angle  ray  portion  of  the  one  hop  F  oblique  echo  trace  is  within  the 
region  of  interest,  the  actual  data  has  a  larger  slope  than  the  region 
of  interest  and  hence  lower  frequency  portions  of  the  high  angle  ray, 
if  they  had  been  recorded,  would  be  excluded  from  further 
processing.  The  retardation  of  the  low  frequency  portion  of  the  low 
angle  ray  of  the  FI  trace  is  caused  by  the  underlying  E  layer.  The 
present  algorithm  will  be  refined  to  include  these  multiple  layer 
effects  at  a  later  time.  Figure  3.6  shows  the  resulting  oblique  F  echo 
trace,  i.e.  only  that  portion  of  the  filtered  oblique  ionogram  data  in 
the  region  of  interest.  Figure  3.7  shows  the  resulting  scaled  oblique  F 
trace  in  white  along  with  the  entire  F  oblique  echo  trace.  In  this  case 
a  simple  algorithm  that  defined  the  scaled  trace  as  the  first  occupied 
range  bin  following  an  unoccupied  range  bin  at  each  sounding 
frequency  was  used. 

The  1800  km  oblique  data  is  seen  in  Figures  3.8  -  3.16. 
The  100  kHz  ionogram,  produced  via  the  frequency  redundancy 
technique  is  seen  in  Figure  3.8.  Figure  3.9  shows  the  resulting  output 
of  the  modified  Wiener  trace  enhancement  and  noise  suppression 
filter.  Figure  3.10  shows  the  regions  of  interest  calculated  from  the  E 
and  F  layers  of  the  corresponding  vertical  ionogram.  The  synthesized 
oblique  F  trace  is  shown  in  blue  in  the  center  of  the  white  region  of 
interest  with  the  synthesized  oblique  E  trace  shown  in  yellow  in  the 
center  of  the  red  region  of  interest.  The  F  region  of  interest  criteria 
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Figure  3.5  Region  of  interest  overlayed  on  the  filtered  oblique  ionogram.  Goose  Bay, 
Labrador  to  Wallops  Island,  Virginia;  October  27,  1988,  1201 UT. 


Figure  3.6  Resulting  oblique  IF  ionogram  in  the  region  of  interest.  Goose  Bay,  Labrador 
to  Wallops  Island,  Virginia;  October  27,  1988,  1201UT. 
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Figure  3.7  Scaled  oblique  IF  trace  shown  in  white  along  with  the  IF  oblique  ionogram. 
Goose  Bay,  Labrador  to  Wallops  Island,  Virginia;  October  27,  1988,  1201 UT. 


Figure  3.8  100  kHz  oblique  ionogram  produced  via  the  frequency  redundancy  technique. 
Goose  Bay,  Labrador  to  Millstone  Hill,  Massachusetts;  April  5,  1988,  1908UT. 
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Figure  3.9  Resulting  100  kHz  oblique  ionogram  after  processing  through  the  modified 
Wiener  trace  enhancement  algorithm  and  the  noise  suppression  filter.  Goose  Bay,  Labrador 
to  Millstone  Hill,  Massachusetts;  April  5,  1988,  1908UT. 


Figure  3.10  Region  of  interest  claculated  from  the  E  and  F  layers  of  the  corresponding 
vertical  ionogram.  Synthesized  oblique  IF  (blue)  and  IE  (yellow)  ionograms  shown  in  the 
corresponding  regions  of  interert.  Goose  Bay,  Labrador  to  Millstone  Hill,  Massachusetts; 
April  5,  1988,  1908UT. 
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Figure  3.11  Region  of  interest  (IF)  overlayed  on  the  filtered  oblique  ionogram.  Goose 
Bay,  Labrador  to  Millstone  Hill,  Massachusetts;  April  5,  1988,  1908UT. 


Figure  3.12  Resulting  oblique  IF  ionogram  in  the  region  of  interest.  Goose  Bay, 
Labrador  to  Millstone  Hill,  Massachusetts;  April  5,  1988,  1908UT. 
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Figure  3.13  Scaled  oblique  IF  trace  shown  in  white  along  with  the  IF  oblique  ionogram. 
Goose  Bay,  Labrador  to  Millstone  Hill,  Massachusetts;  April  5,  1988,  1908UT. 


Figure  3.14  Region  of  interest  (IE)  overlayer  on  the  filtered  oblique  ionogram.  Goose 
Bay,  Labrador  to  Milistone  Hill,  Massachusetts;  April  5,  1988,  1908UT. 
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Figure  3.15  Resulting  oblique  IE  trace  shown  with  the  region  of  interest.  Goose  Bay, 
Labrador  to  Millstone  Hill,  Massachusetts;  April  5,  1988,  1908UT. 


Figure  3.16  Scaled  oblique  IE  trace  shown  in  white  along  with  the  IE  oblique  ionogram. 
Goose  Bay,  Labrador  to  Millstone  Hill,  Massachusetts;  April  5,  1988,  1908UT. 
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are  identical  to  that  used  for  the  2500  km  data.  The  E  trace  region  of 
interest  is  defined  as  all  points  within  three  range  bins  and/or 
sounding  frequencies  of  the  synthesized  oblique  E  trace.  Figure  3.11 
shows  the  F  trace  region  of  interest  overlaid  on  the  filtered  oblique 
ionogram  data.  Note  that  although  the  high  angle  ray  portion  of  the 
one  hop  F  oblique  echo  trace  has  the  same  slope  as  the  corresponding 
portion  of  the  region  of  interest,  the  actual  data  has  an  extended  low 
frequency  segment  of  the  high  angle  ray  that  is  excluded  from 
further  processing.  Figure  3.12  shows  the  resulting  oblique  F  echo 
trace,  i.e.  only  that  portion  of  the  filtered  oblique  ionogram  data  in 
the  region  of  interest.  Figure  3.13  shows  the  resulting  scaled  oblique 
F  trace  in  white  along  with  the  entire  F  oblique  echo  trace.  Figure 
3.14  shows  the  E  trace  region  of  interest  overlaid  on  the  filtered 

oblique  ionogram  data.  Note,  in  this  case  a  segment  of  the  actual  E 
echo  trace,  at  sounding  frequencies  above  the  E  region  of  interest, 

has  already  been  incorrectly  labeled  as  part  of  the  F  trace  and  is 
hence  excluded  from  the  E  trace  at  this  point.  Figure  3.15  shows  the 
resulting  oblique  E  echo  trace,  i.e.  only  that  portion  of  the  filtered 

oblique  ionogram  data  in  the  region  of  interest.  Figure  3.16  shows 

the  resulting  scaled  oblique  E  trace  in  white  along  with  the  entire  E 
oblique  echo  trace. 


4.0  OBLIQUE  IONOGRAM  SCALING  ALGORITHM  STATUS 


The  initial  development  of  the  vertical  end  point 
ionogram-based  oblique  scaling  algorithm  has  been  demonstrated. 
This  algorithm  requires  a  scaled  vertical  ionogram  at  one  end  of  the 
oblique  path  to  generate  an  initial  estimate  of  the  midpoint 
ionosphere  that  is  used  in  the  subsequent  signal  processing  and 
pattern  recognition  stages  of  the  algorithm.  While  this  algorithm 
performs  well  for  relatively  short  oblique  paths,  where  a  reasonable 
estimate  of  the  midpoint  ionosphere  can  be  obtained  from  the  scaled 
vertical  ionogram,  it  performs  in  a  less  reliable  manner  for  cases 
where  the  length  of  the  oblique  path  makes  it  difficult  to  estimate 
the  midpoint  ionosphere  from  the  endpoint  data.  However,  the 
signal  processing  aspects  of  this  algorithm,  which  produce  estimates 
of  the  oblique  echo  traces  from  the  raw  oblique  ionograms,  can  also 
be  utilized  in  the  proposed  model-based  oblique  ionogram  scaling 
algorithm.  A  data  base  of  approximately  200  digital  oblique 
ionograms  has  been  recorded.  These  ionograms  will  also  be  used  in 
the  development  and  evaluation  of  the  model-based  oblique 
ionogram  scaling  algorithm.  The  development  of  the  model-based 
oblique  ionogram  scaling  algorithm,  continued  data  collection  and 
evaluation  of  the  resulting  oblique  ionogram  scaling  algorithm  for  a 
wide  range  of  data  is  in  progress. 


